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The creation of ultra-light dark particles in the late-time FLRW spacetime provides a cosmological 
model in accordance with precise observational tests. The matter creation backreaction implies 
in this context a vacuum energy density scaling linearly with the Hubble parameter H, which is 
consistent with the vacuum expectation value of the QCD condensate in a low-energy expanding 
spacetime. Both the cosmological constant and coincidence problems are alleviated in this scenario. 
We also explore the opposite, high energy limit of the particle creation process. We show that it 
leads to a non-singular primordial universe where an early inflationary era takes place, with natural 
reheating and exit. The generated primordial spectrum is scale invariant and, by supposing that 
inflation lasts for 60 e-folds, we obtain a scalar expectral index n as 0.97. 



I. INTRODUCTION 

The gravitational role of vacuum fluctuations is a chal- 
lenging problem in field theory and cosmology, as vacuum 
is in general not uniquely defined in curved spacetimes 
and its energy density usually depends on the renormal- 
ization procedure. In the case of free (conformal) fields 
in de Sitter spacetime, the renormalized vacuum density 
is A w if 4 jl[ . In a low-energy universe this leads to a too 
tiny cosmological term. In the high energy limit it can 
be used to obtain a non-singular model, with an initial 
quasi-de Sitter universe giving origin to a radiation era, 
through a phase transition in which vacuum decays into 
relativistic particles @- The primordial perturbations 
generated in this process are not scale invariant and a 
subsequent inflation is still needed @. 

In order to see how this non-singular scenario emerges, 
let us take p\ = —A, which is valid in de Sitter space- 
time owing to its maximal symmetry. The conservation 
equation 

p + 3H(p + p)=0 (1) 

assumes the form 

p m + 3H(p m + p m ) = -A, (2) 

where p = p m +A and p = p m — A are the total energy and 
pressure, respectively, and the dot means derivative w.r.t. 
the cosmological time t. We see that a time- varying vac- 
uum density is concomitant with particle production. 

Let us now take A = 3iJ 4 and the relativistic equa- 
tion of state p rn = Pm/i for matter. Eq. (J2J and the 
Friedmann equation 

Pm + A = 3H 2 (3) 

lead to 

H + 2H 2 - 2H 4 = 0. (4) 
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FIG. 1: The Hubble parameter as a function of time for A = 3-ff 4 . 

Apart from the trivial solution H = 1, this equation has 
the solution 

2t = 4 - tanh -1 H, (5) 

where a constant of integration was conveniently chosen. 
This solution, which tends to de Sitter in the asymptotic 
past, is plotted in Fig. 1. We then see that de Sitter 
solution is not stable in the high energy limit, evolving 
to a radiation dominated universe, with vacuum decaying 
into relativistic matter (it is easy to see from §5§ that, 
when H — > 0, we have Ht = 1/2, which characterises a 
radiation phase). As we will see, this process involves 
a backreaction that eventually changes the solution ([3]), 
leading to a subsequent inflation era. 



To clarify these results, let us revisit the seminal work 
by Gibbons and Hawking [|| on the thermodynamics of 
de Sitter spacetime. They show that an observer in a 
de Sitter background experience a thermal bath of tem- 
perature H/2ir, which they have interpreted MS IX ther- 
mal radiation emitted from the observer event horizon. 
The absortion of this radiation by the observer detec- 
tor leads to a shrinking in the horizon area, which is 
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interpreted as resulting from a decrease in the entropy 
beyond the horizon. Nevertheless, that original analy- 
sis does not consider the backreaction associated to the 
particle creation. As discussed above, the corresponding 
energy density and pressure obey, owing to the symme- 
try of de Sitter space, the equation of state p\ = —A, 
where A is the effective cosmological term which defines 
the dc Sitter horizon. Leading this into the conservation 
equation for the total energy ([1]), we obtain equation (|2|). 
As already said, it shows that the particle creation (that 
is, its absortion from the thermal bath) is only possible 
if A decays with time, with a consequent increase in the 
horizon area. This may be understood on the light of a 
reinterpretation of the Gibbons-Hawking thermal bath as 
resulting from the vacuum quantum fluctuations which 
originate the cosmological term. The horizon entropy 
gives the number of degrees of freedom (inside the hori- 
zon) associated to those fluctuations, and its increase just 
means that the extraction of particles from vacuum is an 
irreversible process. This interpretation is in better ac- 
cordance with modern formulations of the holographic 
principle, as we will see below. 

In this paper we shall explore the effect of the parti- 
cle creation backreaction on the solution presented 
above. We will show that it leads to an inflation solu- 
tion with a scale invariant primordial spectrum. Further- 
more, by assuming that inflation lasts for approximately 
60 e-folds, we will be able to obtain the scalar spectral 
index currently observed, n ps 0.97. However, before do- 
ing that, let us consider in the next section the role of 
particle creation in the opposite, low-energy limit of late 
times. We will see that it may also be responsible for the 
present cosmic acceleration, giving origin to a decaying 
vacuum term as well. 



II. THE LATE-TIME LIMIT 

In the case we consider the vacuum energy of inter- 
acting fields, it has been suggested that in a low energy, 
approximately dc Sitter background the vacuum conden- 
sate originated from the QCD phase transition leads to 
A ps m 3 H, where m ps 150 MeV is the energy scale of the 
transition Q . This result is in fact intuitive. In a de Sit- 
ter background the energy per observable degree of free- 
dom is given by the temperature of the horizon, E w H . 
For a masslcss free field this energy is disbributed in a 
volume l/H 3 , leading to a density A ss H 4 , as above. For 
a strongly interacting field in a low energy space-time, on 
the other hand, the occupied volume is 1/m 3 , owing to 
confinement, and the expected density is A ps m 3 H. 

Such a late-time variation law for the vacuum term 
can also be derived as a backreaction of non-relativistic, 
light dark particles production in the expanding space- 
time. Indeed, vacuum fluctuations of a particle of mass 
M < H have a lifetime r ps 1/A/, longer than the Hub- 
ble time, being then observed as real particles irreversibly 
produced at a constant rate r ps 1/r ps M. We can also 



describe this process by saying that, as the wavelength of 
the fluctuations is greater than the horizon radius l/H, 
one of the particles in the pair scapes to the zone be- 
yond the horizon, while the other appears as a created 
particle. If M < H <C 1, the produced particles are non- 
relativistic (photons are not produced in the low-energy 
regime, because there exist no one-loop photon fluctua- 
tions), and the Boltzmann equation for this process is 



where n is the particle number density. By taking p m = 
nM, it can also be written as 

p m + 3Hp m = Tp m . (7) 

The particle production necessarily involves a backre- 
action. Let us take, in addition to ([7]), the Friedmann 
equation 

Pm +A = 3H 2 . (8) 

If our background is approximately de Sitter, the vacuum 
pressure is pa = —A. Using ([7]) and we obtain the 
conservation equation for the total energy, 

P + 3H(p + p)=0, (9) 

provided we take 

A = 2ri? + A , (10) 

where Ao is a constant of integration. Since there is no 
natural scale for this constant, let us make it zero. Then 
we have A = 2TH. This is the time- variation law pre- 
dicted for the vacuum density of the QCD condensate, 
with r ps M ps to 3 . Dividing it by 3 if 2 , we obtain 

v = ^(i-n m )H, (ii) 

where Sl,„ = 1 — fl\ = Pm/(3H 2 ) is the relative mat- 
ter density (for simplicity, we are considering only the 
spatially flat case). In the de Sitter limit {Vt m = 0), we 
have r = 3H/2, that is, the creation rate is equal (apart 
from a numerical factor) to the thermal bath tempera- 
ture predicted by Gibbons and Hawking. It also means 
that the scale of the future de Sitter horizon is deter- 
mined, through T, by the energy scale of the QCD phase 
transition, the last cosmological transition we have. 

For the present time we have from (jlip (with Cl m ps 
1/3) Hq ps r ss m 3 , and hence A ss to 6 , where Hq is 
the current Hubble parameter. The former result is an 
expression of the Eddington-Dirac large number coinci- 
dence [||. The later - also known as Zcldovich's relation 
- gives the correct order of magnitude for A. We also 
have M ps T ps VA ~ 10~ 33 cV. Note that this is the 
quantum of mass expected on the basis of the hologra- 
phy conjecture applied to de Sitter spacetimc. According 
to the holographic principle [H, Q , the number of states 
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inside the dc Sitter horizon is given by the horizon area, 
N pa l/H 2 = 3/ A. On the other hand, the energy in- 
side the horizon is E ps p/H 3 = 3/H. Therefore, the 
energy per degree of freedom is E/N pa vA pa M. As we 
have seen, in the de Sitter limit this quantum of energy 
is equal to the horizont temperature. By the way, let 
us remark that this tiny value for M suggests that dark 
particles may be associated to pseudo-Goldstone bosons 
of the QCD condensate, with a dynamically generated 
mass given by the inverse of the scale of the de Sitter 
horizon. (The cosmological role of an ultra-light pscudo 
Nambu-Goldstone boson with mass ~ Hq was already 
considered before, in another context @). 

The corresponding cosmological model has a simple 
analytical solution, which reduces to the CDM model for 
early times and to a de Sitter universe for t — > oo ■ It 
has the same free parameters of the standard model and 
presents good concordance when tested against type la 
supernovas, baryonic acoustic oscillations, the position of 
the first peak of CMB and the matter power spectrum 
fill [l2l fl3| . Furthermore, as discussed in [l2| , the coinci- 
dence problem is alleviated, because the matter density 
contrast is suppressed in the asymptotic future, owing to 
the matter production. In the right panel of Fig. 2 we 
show the concordance in the Hq x f2 m o space for those 
four observations. We consider the SDSS compilation of 
supernovas la calibrated with the MLCS2k2 fitter, since 
it is less dependent on a fiducial ACDM model than other 
supernovac samples. For the large scale structure (LSS) 
distribution we use the 2dFGRS data. We obtain for 
the present Hubble parameter Ho pa 64 km/s-Mpc, and 
marginalizing over it does not alter the matter density 
concordance value, fi m o ~ 0.45. The best-fit universe age 
is to ~ 14 Gyr. These results do not change if we use for 
supernovas la the Constitution sample compilated with 
MLCS2k2 [ll|]. 

For the sake of comparison we show in the left panel the 
concordance region (with LSS excluded) for the spatially 
flat ACDM model with the SDSS (MLCS2k2) supernovac 
sample. The best-fit Hubble parameter is the same as 
in the A(t) model. However, as already pointed out in 
the literature, in this case there is a tension between the 
SDSS best-fit value for the matter density (ft m .o ~ 0.40) 
and the value obtained from the LSS analysis (Q m o ~ 

0.23) mi. 



III. THE INFLATION SOLUTION 

Let us return to the early-time limit of high energy 
densities. In this case the particles are relativistic, with 
an average energy p pa H. That this is the typical energy 
of any degree of freedom, including vacuum fluctuations, 
can be seen by applying again the holographic principle 
0, H[ to a de Sitter (or quasi-de Sitter) spacetime. The 
number of degrees of freedom inside the dc Sitter horizon 
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FIG. 2: Right. Superposition of the four tests (SNe la, BAO, 
CMB and LSS) performed with the present model. The blue ellipse 
corresponds to the 2<r concordance region. Left. SNe la, BAO and 
CMB 2(7 confidence regions for the spatially flat ACDM model 



is proportional to its surface, N pa l/H 2 , while the total 
energy inside the horizon is E pa p/H 3 pa l/H. There- 
fore, the energy per degree of freedom is p = E/N sa H. 
Another way to support this result is remembering that 
a temperature H/2ir can be associated to the de Sitter 
background. 

The lifetime of a vacuum fluctuation is now given by 
r pa 1 /H . If it is a bit longer than l/H , we have relativis- 
tic particles production at a rate r sa H. The Boltzmann 
equation has now the form 



i|(a*»)=7#». 



(12) 



where n is the particle number density and 7 is a con- 
stant of the order of unity. In an approximately dc Sitter 
spacetime (where p sa constant), (|12[) can be rewritten 
as 



Pn 



3H(p m +p m ) = jHp m , 



(13) 



where p m = p m /3. Equations (|13p and ((7|) are char- 
acteristic of models with interaction in the dark sector 

PI Hi- 

Again, the particle production involves a backreaction. 
From Eq. (fT3"]) (with p\ = —A) and the Friedmann equa- 
tion 



p m + A = 3H 2 , 
we derive the conservation of the total energy, 
p + 3H(p + p) = 0, 

provided we take 

A = ^ 2 + A . 



(14) 
(15) 

(16) 



Since there is no natural scale for the integration con- 
stant A , we take it equals to zero. Let us remark that a 
quadratic dependence of A on H was already considered 
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in other contexts. For example, it can be obtained on the 
basis of renormalization group results (l6| . 

The solution for ff!3]) . (fT4"|) and (HU) is given by 



1 

H= (f 



1 



(17) 



where e — (4 — 7)/2 and _ff is now an integration con- 
stant. For 7 w 4 we have e <C 1. Hence, for e\t\ <C 1 
we have H ss ifoj that is, an approximately dc Sitter 
solution, as assumed above. 

We can verify from (fT?]) that e is the slow-roll inflation 
parameter defined by e = —H/H 2 . As it is constant, 
the second slow-roll parameter is 5 = — e, and the scalar 
spectral index is [13] 



n. = 1 — 2e. 



(18) 



From (|17j) we define the time when inflation ends by 
Hot = 1/e. Therefore, assuming that inflation lasts for 
about 60 e-folds, as required to solve the horizon prob- 
lem, we have e « 1/60, leading to n s sa 0.97. This is 
approximately the current value obtained for the scalar 
spectral index from CMB observations, n s = 0.96 ± 0.01 
(la) [l8|]. Note that such a relation between the number 
of e-folds and n s is not general, holding only for models 
with constant e. On the other hand, it is easy to show 
from (Q3J) and (fll)|) that the matter relative density is 
&m — e/2 ~ 0.01. That is, during inflation radiation and 
matter represented around 1% of the energy content. 

Up to now we have been considering 7 as a constant 
in (fT2|) . However, it actually decreases as long as H 
decreases, because the production of a particle species 
stops when H crosses the threshold of the particle mass 1 . 
Eventually 7 goes to zero, hence we have e = 2 and, for 
et > 1/H , it follows from ([17]) that Ht = 1/2, which 
characterises a radiation phase, with O m = 1. As for the 
ultra-light dark particles considered above, their produc- 
tion starts after the QCD phase transition, but they are 
subdominant during the radiation phase [101 ] . 

The energy scale of inflation (given by the integration 
constant Ho in (fTT|) ) is not determined in this model. It 
is reasonable to assume that it is above the electroweak 
scale, in which case all the observed particles are pro- 
duced. The Planck scale (H « 1) is a natural superior 
limit, at which the ansatz A ss H discussed at the be- 
ginning is valid. In this case inflation would take place 
at the end of the non-singular quasi-de Sitter phase. 

As in this model we have two interacting fluids, namely 
the vacuum term and relativistic matter, the presence 
of entropic perturbations must be investigated. Suppos- 
ing that both components have no intrinsic non-adiabatic 



1 Strictly speaking, the temperature and entropy associated to the 
horizon are not well defined when the spacetime departures sig- 
nificantly from do Sitter. 




per turbations, the only entropic contribution is given by 

m 



(19) 



where a hat means perturbation of a given quantity. 

Owing to its equation of state, p\ = —A, the vacuum 
component perturbation is negligible for scales inside the 
horizon [Hjj . Therefore, the relative entropic perturba- 
tion (IT9l reduces to 



p A 

P~-P~ -pr> 
P P 



Pi, 



(20) 



where in the last approximation we have used the fact 
that, during the inflationary phase, A f=a p. The non- 
adiabatic perturbations are then proportional to p m , and 
their presence does not affect the scale-invariance of the 
primordial espectrum. 

In reference 0], both the Planck limit Awff 4 and the 
low-energy result A w MH were also modeled, alterna- 
tively, by a self-interacting scalar field, with the interac- 
tion potential V(cf)) playing the role of a decaying vacuum 
term, and the kinetic term <fi 2 /2 representing produced 
stiff matter. Let us show that we can do the same here, 
that is, we can model this inflationary scenario by intro- 
ducing an equivalent inflaton field. The field equations 
are Q 



3H 2 



V + 2H' 2 , 



= -2H' 



(21) 



(22) 



where the prime means derivative with respect to 
From equation (|T6|) we see that, in this case, V 
3 r yH 2 /4. The solution is then given by 



1 

H 



v(4>) = 



et 



3lH 2 



1 

,±V2e<t 



(23) 



(24) 



where now e = (3/2)(4 — 7)/2. The first equation is, 
apart a re-scaling of time, the same solution (fT7| obtained 
before. The second equation gives, for e <C 1 (that is, 
7 w 4), an approximately constant potential V ~ 3H 2 . 
It is easy to verify that V"/V = (V /V) 2 = 2e < 1. 
This guarantees that we have a slow-roll inflation, with 
a scale-invariant spectrum k 3 P = H 2 /9e [l?} • It also 
means that, as asserted above, the perturbations in the 
potential (i.e., in the vacuum term) are negligible, since 
V = iH 2 ^r± < 4>. The potential JUJ was also con- 
sidered by [l9j in the context of a model with particle 
creation. 
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IV. CONCLUSION 

Particle creation is something expected in expanding 
spacetimes (2Cj . In spite of the difficulty in deriving the 
production rate and backreaction in general, this phe- 
nomenon may in principle be related with inflation and 
the present cosmic acceleration, a possibility already con- 
sidered in different ways by some authors [U, [22j . We 
have shown in this paper that the production of non- 
relativistic, ultra-light dark particles at late times leads 
to a concordance model in which the vacuum density de- 
cays with time. In the limit of high energies, on the other 



hand, we have shown that the production of relativistic 
particles from the vacuum leads to a viable inflationary 
solution. 
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